A phonon laser in ultra-cold matter 
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We show the possible excitation of a phonon laser instability in an ultra-cold atomic gas confined 
in a magneto-optical trap. Such an effect results from a negative Landau damping of the collective 
density perturbations in the gas, leading to the coherent emission of phonons. This laser instability 
can be driven by a blue-detuned laser superimposed to the usual red-detuning laser beams which 
usually provide the cooling mechanism. Threshold conditions, instability growth rates and satura- 
tion levels are derived. This work generalizes, on theoretical grounds, the recent results obtained 
with single ion phonon laser, to an ultra-cold atomic gas, where real phonons can be excited. Future 
phonon lasers could thus adequately be called phasers. 

PACS numbers: 



The last decade the concept of phonon laser has been 
studied in several different systems, e.g. quantum wells 
[I], ions 0, nanomechanics 4 , and nanomagnets [S]. 
Furthermore, very recently the possibility to create a 
Doppler instability and create a phonon laser with a 
single trapped atom or ion has been considered [5J [7j. 
Here we show that it is possible to generalized the single 
trapped atom configuration to the case of an large ensem- 
ble of atoms, thus enabling collective phonon excitations. 
In particular, the excitation of a phonon laser instability 
in an ultra-cold atomic gas confined in a magneto-optical 
trap (MOT) is described. Such an effect results from a 
negative Landau damping of the collective density per- 
turbations in the gas, leading to the coherent emission of 
phonons. In contrast with the single ion case, here the 
phonon frequency is not determined by the single atom 
oscillation frequency due to the parabolic confining po- 
tential, but by the boundary conditions of the internal os- 
cillations of the gas. It therefore corresponds to a collec- 
tive oscillation, instead of a single particle effect. In that 
respect, it is closer to the paradigm of an optical laser, 
where the photon modes correspond to internal vibra- 
tions of the optical cavity, thus making the name phaser 
a natural choice for such collective modes. The acous- 
tic oscillations can also, in principle, be coupled to the 
outside world, by mechanical or electromagnetic means, 
thus providing a source of coherent acoustic radiation. 

Collective processes in ultra-cold gas clouds have been 
considered by several authors in recent years. These pro- 
cesses result from the existence of an effective atomic 
charge in the gas cloud 8 , which create collective forces 
of the Coulomb type [9] . It was also shown that plasma- 
acoustic waves can be excited in the gas, which are similar 
to the usual sound waves but with a lower frequency cut- 
off [TU] ■ These waves can provide the physical support of 
our laser instability. 

Like in the case of a single trapped ion, the collective 
laser instability can be driven by a blue-detuned laser, su- 
perimposed to the usual red-detuned laser beams which 
usually provide the cooling mechanism. But in contrast 



with the single ion phonon laser, real phonons are excited 
in the medium and can eventually be coupled to an exte- 
rior system, just like photons produced in current optical 
laser cavities. 

The dynamical behavior of the ultra-cold atomic gas 
will be described by a wave kinetic equation, where the 
atom recoil effects due to the emission and absorption 
of both photons and plasmons are retained llj. This 
equation corresponds to a straightforward generalization 
of well known quantum kinetic equations |12| 1X3] , which 
include a collective (or mean-field) potential. The insta- 
bility growth rates are derived, in the linear regime. Sat- 
uration of the phonon laser instability is also described 
with the help of a quasi-linear version of the wave kinetic 
equation. Threshold conditions for the occurrence of the 
phonon laser instability are also established. This work 
therefore extends a recently published concept of single 
ion phonon laser [7] into the collective regime, by con- 
sidering an ensemble of ultra-cold atoms confined in a 
MOT. In contrast to previous results, real phonons can 
be excited and eventually coupled to external environ- 
ments. 

The starting point for our analysis is provided by the 
quantum wave kinetic equation for the ultra-cold gas 
jTTJ [T2] • This equation describes the space and time 
evolution of the Wigner function W(r, v, t) associated 
with the kinetic state of the atoms in the cloud, where 
r and v are the position and the velocity of the atomic 
centre-of-mass. It is well known that this quantity pro- 
vides the quantum analogue to the classical representa- 
tion of the atom centre-of-mass in phase space. It is 
normally referred to as a quasi-probability distribution 
since it can take negative values. Here we retain the ex- 
act quantum description of the centre-of-mass state of 
motion, which means that all the recoil effects associ- 
ated with the atom interaction with electromagnetic and 
acoustic fields are taken into account. We will see that 
recoil is an essential ingredient of our phonon laser ef- 
fect, but the classical limit will also be discussed. The 
evolution equation for the quantum quasi-probability dis- 
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where the total time derivative is <i/<ii = (d/dt + v • V). 
This equation states that the changes in the distribution 
W(r,v,t), are due to three different factors. The first 
is the external force term, associated with the collective 
mean-field potential as defined by 
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where W(k) {±) = W(v ± Kk/2M), with M the mass of 
the atom, and we have used a spatial Fourier decom- 
position in the integral in This collective force de- 
scribes the atom-atom interactions inside the gas and is 
determined by the potential V, formally analogous to an 
electrostatic potential, governed by the Poisson equation 
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Here Q determines the effective charge of the atom 
[HI |9j [10], as will be discussed later, and depends on 
the intensity of the incident cooling and pumping laser 
beams. The second source in Eq. (JTJ is due to the red- 
detuned cooling laser beams (normally six in number) 
with electric fields of the form E(r, t) = E„ exp(zk„ • r — 
iu> a t), corresponding to the same frequency u> a , which is 
nearly equal but slightly lower than the frequency of the 
atomic transition used for laser cooling, and wavenum- 
bers |k„| = k a = uj a /c. In equation 0, the cou- 
pling coefficients g n are associated with the red-detuned 
laser beams. Finally the third source term in Eq. (JT|), 
is due to the blue-detuned laser pump field, E(r,i) = 
E(,exp(ikb • r — ioj^t), with frequency cjf, > tu a . This 
pump field will produce a population inversion in the 
center of mass velocity state, leading to an instability 
of the acoustic waves in the ultra-cold gas. We already 
now stress that these acoustic waves have some special 
features which make them clearly distinct from common 
acoustic oscillations. This will become apparent later. 

It is obvious from equation (JlJ that the first term cou- 
ples the atomic states with the spectrum of possible low 
frequency oscillations of the collective potential ^(k), 
with |k| <§C kf, | ~ |k„|, and the second term couples the 
atoms with the blue-detuned pumping field. These two 
terms can describe the atomic recoil associated with the 
exchange of momentum with both photons (second term) 
and phonons (first term). The radiation coupling coeffi- 
cient, which is the mediator of photon recoil, is given by 
the usual expression g = — (E / fij^^i^i) , where d 2 i is 
the dipole matrix element of the radiative transition be- 
tween the internal atomic states 1 and 2, and P21 is the 



density matrix element, assumed nearly at equilibrium 
with the radiation field. In the above kinetic equation 
we have neglected the confining potential, which is as- 
sumed to be parabolic. This is justified by the fact that 
the wavelength of the collective oscillations are assumed 
much shorter than the size of the confining potential well. 
Therefore, as a first approximation, the medium can be 
assumed as infinite. The role of the boundary conditions 
will be discussed at the end. 

In order to simplify the above description we consider 
an one-dimensional model, by assuming v — u(k/k)-\-~v±, 
and integrate over the perpendicular velocities with re- 
spect to the pump laser beam. In the low frequency spec- 
trum of the collective potential l^(k), we retain spectral 
components propagating in the same direction, k || kf,. 
We define the quasi-distribution for the parallel veloci- 
ties as G(u) = J W(u, v±)dv±. The resulting evolution 
equation is 
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where v is a phcnomcnological viscosity coefficient. In 
reality, this quantity depends on the angle between the 
pumping and the cooling beams. But, for order of magni- 
tude estimates, we can use the the viscosity coefficient, as 
defined in the absence of the pumping laser beam [9j [10] . 
The distribution Go(u) is that of the laser cooled gas, 
as obtained in the absence of the pumping laser. This 
last term in Q plays the role of spontaneous emission 
by depopulating the higher velocity states driven by g, 
as shown below. 

Three physical processes are assumed to occur in par- 
allel, as described in figure 1. We can see a very strong 
resemblance with the usual three level laser model. First, 
the red-detuned laser beams, will create a very low tem- 
perature quasi-distribution G(u), which corresponds to 
our ground state. Second, the blue-detuned pump laser 
field will excite high velocity atomic states around the 
velocity u + hk^/M, due to the absorption of photon 
momentum, therefore populating the upper level for the 
phonon laser transition. This high velocity distribution 
corresponds to a kind of particle beam, thus creating a 
population inversion in the centre-of-mass states. It will 
then lead to a negative Landau damping of acoustic- 
like oscillations with frequency to and wavenumber fc, 
which will result in the coherent emission of phonons. 
Due to phonon emission, a third and intermediate ve- 
locity state will be populated. The temporal evolu- 
tion of the phonon field will be dictated by a com- 
plex frequency u> = uj r + i{^k — v 1^)- We assume that 
G(u) = G(u) + G(u), where G(u) is the equilibrium dis- 
tribution, and G(u) is the perturbation, which is assumed 
to evolve as exp(ik • kr — icot). By linearizing equations 
([lj and ([3]), it is possible to show that the instability 
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FIG. 1: Phonon laser scheme: (a) The gas is cooled down by 
red-detuned laser beams with frequency u a and wavevector 
k a into a velocity distribution Wo(v); (b) the ultra-cold gas 
is pumped by a blue-detuned laser with frequency ui and 
wavevector kt, , therefore creating an inversion of population in 
velocity space; (c) Phonons with frequency u) and wavevector 
k are coherently emitted by the atoms, which decay into a 
lower kinetic energy state. A one-dimensional configuration 
is assumed for simplicity. 



growth rates are determined by [TT] 
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where the population difference is calculated for the res- 
onant parallel velocity u r = k juj r , and the acoustic oscil- 
lations obey the dispersion relation 

(6) 
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This dispersion relation describes a modified form of 
acoustic oscillations. The first term on the right hand 
side of Eq. (|6| corresponds to a plasma frequency cut- 
off, the second term is the usual acoustic dispersion term, 
and the last term is a purely quantum correction due to 
the atom recoil. The plasma frequency cu p and the ion 
sound speed u s are determined by 
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where n is the unperturbed gas density. By comparing 
this definition of the plasma frequency of the neutral gas 
with the usual definition valid for an electron-ion plasma, 
we conclude that the effective charge of the neutral atoms 
inside the ultra-cold gas is y/eoQ, as first considered by 
[5]. It is clear that we will have an acoustic wave growth 
if the inverse Landau damping is positive 7^ > 0, and if 
it is large enough to compensate for the wave losses, i.e. 
7^ > v. We can then write the threshold condition as 
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FIG. 2: Saturation of the laser instability, as described 
by the quasi-linear equations. Temporal evolution of the 
normalized population inversion in velocity space z(r) = 
AG(u, t)/AG(u, 0), and the phonon amplitude square y(r) = 
I Vfc (t) 1 2 , as a function of the time variable r, for a = 1 and 
v — 1/2. For illustration purposes, we take y(0) = 0.001 



It should be noticed here that such an instability remains 
in the classical limit. In this limit, we can neglect the 
atom recoil due to emission or absorption of phonons, 
which allows us to develop G^\ around G(u), as G^ ~ 
G(u) ± (hk/M)(dG/du). The threshold condition now 
reads 



dG\ 2uMk 2 
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The population inversion is now represented by the 
derivative of the atom distribution at the resonant ve- 
locity u r — uj r /k. The resemblance with a three level 
laser will be somewhat lost, but the physical principle 
stays the same. Next, we discuss the instability satu- 
ration. The acoustic wave growth will tend to decrease 
the population difference between the upper and lower 
kinetic energy levels. This process is approximately de- 
scribed by the quasi-linear equations |llj 



dG _ |y(k,t)| 2 
dt ~ W- 



[G(u - hk/M) - G(u + hk/M)] (10) 



In order to be specific, we assume that the mean distribu- 
tion G{u) is divided into three distinct regions of spectral 
interest, centered arround uq — 0, u\ — h(kt — k)/M and 
u 2 = hk b /M. We have, G(u) = G (u) + Gi(u) + G 2 (u). 
For the extreme case of Gj(u) = GjqS(u — Uj), with 
j = 0, 1,2, the above quasi-linear equation would be re- 
duced to simple balance equations between the three ve- 
locity level populations. The resonant atom velocity de- 
fined in equation (|5| for the inverse Landau damping will 
then be exactly equidistant from U\ and u 2 - Integrating 
over the parallel velocity, and defining the population dif- 
ference SG — jj\G 2 (u) — G\(u)\du, for the atoms centered 
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around the upper and intermediate parallel velocities u\ 
and U2, we can easily establish the evolution equation 
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On the other hand, the inverse Landau damping of the 
phonon mode will also change and, as a consequence, the 
phonon square amplitude (energy) will evolve according 
to 
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where 7^ is the growth rate defined above. Introducing 
a new time variable, r, and a new viscosity coefficient v, 
such that 

r = (hk 3 /7rQcj r )AG(0)t , v = {r/t)v (13) 
we can rewrite these two equations as 
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with z ee AG/AG(0), y = 
(irQui r /fi 3 k 3 ). These coupled evolution equations clearly 
exhibit the threshold condition z(0) > 2D, already stated 
in equation (|8j|, and saturation for increasing values of 
time. This is illustrated in figure 2. We can see that the 
population difference decreases along time, due to the in- 
crease of the phonon coupling between the two velocity 
states around u\ and u-i. As a consequence, the growth 
rate slows down, leading asymptotically to saturation. 

In our discussion we have only considered a single 
phonon mode, corresponding to a given wavenumber k. 
It is important to discuss the selection mechanisms iden- 
tifying such a mode inside the phonon spectrum. Several 
different processes can lead to mode selection. First, we 
have spontaneous selection of the mode with the largest 
growth rate. Due to the exponential growth, such a mode 
can easily detach from the background noise. However, a 
more effective mechanism is due to the finite size of the 
gas cloud, which can be confined in a MOT. This leads 
to a discretization of the phonon mode spectrum, as re- 
cently discussed [10j . and the resulting discrete phonon 
modes with a well defined internal structure are called the 
Tonks-Dattner resonances of the nearly spherical gas. In 
a certain sense, this can be seen as the phonon equivalent 
of the laser modes in a typical laser cavity. The domi- 
nant phonon mode in an ultra-cold gas will be due to the 
Tonks-Dattner resonance mode with the largest inverse 
Landau damping. Recent experimental work suggests 
that such modes can be spontaneously excited in the ab- 
sence of any pump laser beam |14j . Typical experimental 



conditions are 10 9 — 10 10 confined atoms, and plasma fre- 
quency (the lower cut-off for Tonks-Dattner oscillations) 
of 100 to 200 Hz. Finally, these modes can be excited by 
an external source, for example due to amplitude modu- 
lation of one of the incident laser beams, and the phonon 
laser system will operate as an amplifier. Such a phonon 
laser system, operating as an oscillator or as an amplifier, 
could adequately be called a phaser (instead of saser, as 
proposed by others). 

In conclusion, we have shown that a phonon laser in- 
stability can be excited in the ultra-cold atomic cloud 
confined in a MOT, and pumped by a blue-detuned laser 
beam. Threshold conditions, linear growth rates and 
nonlinear saturation have been established. We have 
used a quantum kinetic description, which retains the 
atom recoil effects associated with the emission and ab- 
sorption of both photons and phonons. This work gen- 
eralizes recent proposals for a single atom phonon laser 
[IH |7] to a large ensemble of atoms . The present phonon 
laser configuration is able to produce real phonons in the 
gas. Due to the acoustic cavity, the phonon modes are 
determined by the boundary conditions of the gas, in 
contrast to the single atom oscillation frequency studied 
in previous works. 
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